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VANADATE-SULFATE  MELT  THERMOCHEMISTRY  RELATING  TO 
HOT  CORROSION  OF  THERMAL  BARRIER  COATINGS 


INTRODUCTION 

The  purpose  of  this  Report  is  to  draw  together  the  knowledge  of  vanadate-sulfate  thermo¬ 
chemistry  relating  to  TBC  hot  corrosion,  including  that  obtained  in  work  done  over  several  years  in 
Code  6170  of  the  U.S.  Naval  Research  Laboratory,  and  to  present  the  total  in  reviewed  form.  This  is 
important  for  several  reasons.  First,  gas  turbine  technology  is  now  moving  toward  the  wide-spread 
use  of  TBCs  on  the  1st  stage  turbine  vanes/blades  (i.e.,  highly  critical  components),  with  many  such 
TBC-fitted  engines  likely  to  be  used  for  electricity-generation  or  other  industrial/marine  applications 
where  V,S-containing  fuels  are  possible.  The  hot  corrosion  of  TBCs,  and  therefore  the  thermo¬ 
chemistry  of  vanadate-sulfate  melts,  may  become  of  concern  in  these  cases.  Second,  although  there 
have  been  many  thermodynamic  studies  of  metallurgical  slags,  fused  salts,  silicate  glasses,  etc., 
surprisingly  few  investigations  of  vanadate-sulfate  melts  have  been  made,  despite  the  long  history  of 
vanadate-sulfate  hot  corrosion  (1).  And  third,  among  even  these  few  studies,  significant  differences 
in  findings  have  been  reported.  By  defining  these  differences,  it  may  be  possible  to  identify  areas  In 
vanadate-sulfate  thermochemistry  where  new  research  is  needed. 


Gas  Turbine  Vane/Blade  TBCs 


Zirconia  is  a  high  temperature  thermal 
insulator,  and  a  ZrOj  "thermal  barrier  coating"  as 
thin  as  0.010"  (0.25  mm)  can  reduce  the  temp¬ 
erature  between  the  engine  gas  and  component 
metal  by  1 0O-ZOO’C.  Many  problems  have  had 
to  be  solved  to  make  TBCs  viable  (2,3),  espec¬ 
ially  for  1  St  stage  turbine  vanes  (stationary)  and 
blades  (rotating)  where  large-scale  TBC  failure 
could  jeopardize  engine  life.  The  TBC  structure 
now  commonly  used  for  vane/blade  applications 
(Fig.  1 )  consists  of  an  0.003-0.005"  oxidation- 
resistant  metallic  "bond  coat"  of  MCrAlY  (where 
M  =  Co,  Ni,  or  NiCo)  or  Pt-aluminide  under  a 
-0.010"  layer  of  Y203(8-wt%)-Zr02  having  a 
columnar,  carpet  pile-like  physical  structure. 
This  columnar  structure  is  produced  by  electron 
beam-physical  vapor  deposition  (EB-PVD),  and  Is 
thought  to  be  essential  for  relieving  thermal 
stress  within  the  TBC,  and  so  yielding  long 


Fig.  1  Schematic  drawing  of  EB-PVD  thermal 
Significant  advantages  are  claimed  for  barrier  coating. 

TBCs  on  vanes/blades,  including  improved  power 

(up  to  20%),  fuel  economy  (several  percent),  and/or  component  life  (up  to  3X  improvement).  For 
example,  in  land-based  gas  turbines,  an  increase  of  55*C  in  turbine  gas  temperature  Is  expected  to 
provide  an  8  to  1 3  percent  gain  in  power,  and  a  1  to  4  percent  increase  in  simple  cycle  efficiency  (4). 
The  faith  of  the  gas  turbine  Industry  In  TBC  technology  is  shown  by  the  fact  that  rnany_$JIOg_gf 
millions  of  dollars  have  been  spent  by  such  companies  as  United  Technologies  (Pratt  &  Whitney 
Aircraft  Engines),  General  Electric  Aircraft  Engine  Company,  Chromalloy  Corp.,  and  Praxair  Surface 
Technologies  to  purchase  production-scale,  EB-PVD  coaters  for  application  of  TBCs  to  aviation  and 
industrial  gas  turbine  vanes/blades. 
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Combined  Cycle  Gas  Turbines  for  Electricity  Generation 


A  world-wide  movement  is  underway,  as  exemplified  in  the  U.S.  Dept,  of  Energy  8-year, 
$700M  Advanced  Turbine  Systems  (ATS)  Program  (5),  to  use  high  efficiency  gas  turbines  in  combined 
cycle  with  bottoming  steam  turbines  to  produce  electricity  with  up  to  60%  thermal  efficiency.  Such 
plants  would  require  approximately  50%  less  fuel  than  present-day  fired  boiler  plants  which  have  35- 
42%  thermal  efficiency.  Moreover,  since  these  plants  are  Intended  mostly  to  burn  natural  gas,  which 
normally  consists  of  90%  or  more  of  methane  (CH4),  there  would  be  proportionately  less  CO2  released 
to  the  environment  than  with  fuels  of  higher  C/H  ratio. 


To  achieve  high  fuel  efficiency,  the  massive  electricity-generating  gas  turbines  (with  power  In 
ranges  of  up  to  250-300  MW)  will  employ  the  latest  gas  turbine  technology,  with  many  using  single 
crystal  blades  (despite  their  large  size)  and  TBC  coated  vanes/blades.  Because  of  the  possibility  of 
future  high  prices,  or  shortages,  in  natural  gas,  these  engines  must  have  the  back-up  capability  for 
firing  with  industrial-grade  petroleum  fuel.  The  vane/blade  TBCs  must  therefore  have  the  ability  to 
withstand  at  least  certain  levels  of  vanadium  and  sulfur  fuel  impurities,  even  though  the  high  surface 
temperatures  for  TBCs  may  tend  to  reduce  deposit  accumulation.  Also,  it  is  common  in  less  developed 
countries  to  burn  locally-produced  industrial  or  crude  oils  that  may  contain  significant  vanadium  or 
sulfur,  and  engine  manufacturers  are  seeking  vane/blade  TBCs  that  can  be  used,  with  economic 
advantage,  even  in  these  circumstances. 

Vanadium  in  Petroleum 

Vanadium  occurs  in  most  crude  oils,  and  in  some  cases,  at  levels  up  to  500  ppm  (6).  The 
vanadium  is  mostly  complexed  as  heavy  metalloporphyrin  molecules,  and  tends  to  remain  in  the 
residual  during  refining;  however,  since  metalloporphyrins  are  slightly  volatile  at  higher  temperatures, 
some  vanadium  may  also  be  found  in  the  higher-boiling  distillate  fuels. 

Vanadium  (as  VjOg)  accumulates  on,  and  deactivates,  the  catalysts  used  in  the  fluid  catalyst 
cracking  of  petroleum.  Much  effort  has  been  spent  to  alleviate  this  problem,  which  is  now  mostly 
approached  by  catalyst  modification,  or  metals  passivation  where  certain  metals  such  as  Sb  or  Sn  are 
added  to  the  feedstock  to  "passivate"  the  catalyst  against  V^Oe  attack  (6,7).  Vanadium  also  interferes 
with  the  hydrogen-desulfurization  (HDS)  of  high-sulfur  crudes  by  depositing  on  the  catalyst.  To  remove 
vanadium  Itself,  refiners  can  use  an  aggressive  hydrogen  treatment  over  catalysts  (hydrodemetallization 
or  HDM),  as  well  as  coking  or  other  means  (6).  Approximately  2,650  metric  tons  of  vanadium  were 
recovered  and  marketed  in  the  US  from  petroleum  residues  in  1996  (8).  The  treatments  for  removing 
vanadium  are  sophisticated  and  expensive,  however,  and  may  not  be  routinely  performed  in  refining, 
especially  in  less  developed  countries. 

Previous  Studies  of  Vanadate  Hot  Corrosion  of  TBCs 

the  fuel  crises  of  the  1 970s,  the  U.S.  Dept,  of  Energy  investigated  alternative  sources 
of  liquid  fossil  fuel  (coal  liquefaction,  etc.),  as  well  as  means  for  improving  engine  fuel  economy  (e.g., 
thermal  barrier  coatings).  The  possible  corrosive  attack  of  TBCs  by  low-quality  alternate  fuel  con¬ 
taining  vanadium,  phosphorus,  sulfur,  etc.  was  investigated  under  DOE  contract  by  several  engine 
companies  including  Westinghouse,  General  Electric,  Solar  Turbines  and  others.  The  results  from  these 
studies  were  primarily  published  in  two  Conference  Proceedings  (9,10)  jointly  sponsored  by  DOE  and 
EPRI  (Electric  Power  Research  Institute). 

The  DOE-sponsored  studies  focused  essentially  on  Mg0(20-25wt%)-Zr02  (MSZ)  and  Y^OafS- 
20wt%)-ZrO2  (YSZ)  thermal  barrier  coatings.  MSZ  was  found  to  lack  high  temperature  phase  stability, 
and  to  be  destabilized  simply  by  thermal  cycling  to  1000”C.  The  MgO  stabilizer,  although  somewhat 
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more  resistant  to  reaction  with  VjOg  than  Y2O3,  reacted  readily  with  SO,  to  form  MgS04,  ^ftd  was 
strongly  leached  from  MSZ  by  the  SO,  in  burner  rigs  burning  1  wt%-S  fuel.  YSZ  was  clearly  superior, 
and  exhibited  good  1 000*0  phase  stability,  as  well  as  resistance  to  reaction  with  Na2S04-S0,  at  the 
SO,  partial  pressures  found  in  engines  burning  fuel  of  1-2%  S  levels.  However,  Y2O,  was  confirmed 
by  many  as  being  highly  reactive  with  VjOg.  even  when  the  Y2O3  was  at  reduced  activity  as  8-wt% 
(4.5  mol%)  Y2O3  In  solid  solution  in  ZrO,. 

Building  upon  the  original  DOE-funded  work,  the  U.S.  Naval  Sea  Systems  Command  sponsored 
research  to  explore  the  potential  of  TBCs  for  vane/blade  use  in  the  General  Electric  LIVI2500  gas  turbine 
which  powers  several  important  classes  of  Navy  ships.  Burner  rig  tests,  along  with  one  limited  engine 
test.  Indicated  that  coatings  of  Y203(8-wt%)-Zr02  in  fact  resisted  NaS04-S03  hot  corrosion  better  than 
the  metallic  MCrAlY  or  Pt-alumlnide  coatings  normally  used  for  corrosion  protection  (11).  There  was 
also  no  evidence  of  molten  salt  penetration  into  the  intergranular  spacings  of  the  EB-PVD  prepared 
thermal  barrier  coatings.  The  relative  performance  of  MgO-,  CeO,*,  and  Y203-stabillzed  ZrOj  TBCs  was 
evaluated  in  another  burner  rig  test  series,  which  used  fuel  contaminated  with  vanadium  (up  to  ~90 
ppm),  sea  salt,  and  sulfur  (12).  MgO-ZrO,  was  judged  to  be  by  far  the  most  reactive,  whereas  CeOj* 
ZrO,  and  Y203-Zr02  were  judged  to  perform  moderately  well  (in  terms  of  the  TBC  not  spalling).  All 
showed  chemical  reactivity,  however,  with  MgO  forming  large  amounts  of  MgS04,  and  CeO,  and  Y2O3 
each  forming  vanadates,  although  with  some  sulfate.  The  formation  of  cerium  vanadate  (CeV04)  was 
somewhat  unexpected,  since  CeO,  shows  some  resistance  to  V2O6  reaction  (see  below),  but  this  may 
have  resulted  because  of  the  relatively  high  concentrations  of  vanadium  and  sulfur  used. 

THERMOCHEMISTRY  OF  VANADATE-SULFATE  MELTS:  REVIEW  AND  DISCUSSION 

The  principal  driving  force  for  the  degradation  of  TBCs  by  hot  corrosion  is  chemical  reaction 
between  the  stabilizing  oxides  (e.g.,  Y2O3,  CeO,/  etc.)  and  the  molten  vanadate-sulfate  engine  deposits, 
or  essentially  vanadate-sulfate  melts.  Vanadium-containing  melts  of  relatively  low  V2O6  activity  (e.g., 
NaVO,)  can  cause  phase  reversion  in  Ce02-Zr02  TBCs  (and  possibly  In203-Zr02  TBCs)  by  a 
“mineralization"  effect  where  no  chemical  reaction  of  the  CeO,  can  be  detected  (13).  However,  this 
is  of  second  magnitude  compared  to  hot  corrosion  chemical  reaction.  The  thermochemistry  of  the 
vanadate-sulfate  melt  determines  whether  chemical  reaction,  or  corrosion,  of  the  TBC  components  will 
occur,  and  it  Is  thus  the  ultimate  factor  in  deciding  TBC  life  in  corrosive  environments. 

Mechanism  of  Vanadate  Hot  Corrosion  of  Y.>0,-Zr0,  TBCs 

As  established  in  DOE-sponsored  research  (9),  and  since  verified  by  others,  the  attack  of 
vanadium  on  yttria-stabilized  zirconia  (YSZ)  thermal  barrier  coatings  occurs  by: 

YjO,  (in  solid  soln  In  ZrO,)  +  V2O6  (D  =  2  YVO4  (s)  +  ZrO,  (s,  monoclinic)  (11 

where  Y2O3  in  solid  solution  In  the  ZrO,  matrix  reacts  with  the  V2O6  (I)  component  of  the  vanadate- 
sulfate  engine  deposit  to  produce  highly  stable,  solid  YVO4  (mp  1810*0  and  solid  ZrO,,  but  with  the 
ZrO,  being  destabilized  and  in  the  monoclinic  phase  structure.  No  chemical  reaction  of  VjOg  with  ZrO, 
is  normally  detected.  Reaction  (1 1  can  be  confirmed  by  heating  thin  deposits  of  NaVO,  on  sintered  YSZ 
pellets,  where  one  then  sees  masses  of  acicular  YVO4  crystals  formed  on  the  YSZ  surface  by  outward 
diffusion  of  YjO,,  while  destabilized  monoclinic  ZrO,  is  left  below  (14).  Specific  details  of  the 
crystallographic  orientation  of  YVO4  growing  on  YSZ  have  been  revealed  by  thin  film  TEM  (1 5),  while 
the  diffusion  of  VjOg  and  YjO,  within  YSZ  has  been  studied  by  Rutherford  back-scattering  (16). 

The  activity  of  VjOg  (I)  determines  whether  reaction  (11  will  occur,  since  YVO4  and  ZrO,  are 
pure  solids  (activity  =  1),  and  although  Y203is  in  solid  solution  In  ZrO,  (act.  <  1),  its  activity  is  fixed 
until  reaction  (1 1  actually  commences.  It  is  important  therefore  to  know  the  activity  in  vanadate- 
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sulfate  engine  deposits  to  predict  when  reaction  [11  is  possible,  and  to  evaluate  alternative  stabilizers 
for  improving  the  hot  corrosion  resistance  of  TBCs.  As  an  example,  the  activity  coefficient  (y)  for 
0.045  mol-fraction  (8-wt%)  of  YjOj  in  Zr02  ZSOO’C  has  been  measured  as  ~0.1  (17),  and  is 
probably  lower  at  BOO-SOO'C  (assume  0.01),  since  negative  deviation  from  ideal  solution  behavior 
normally  increases  with  decreasing  temperature.  Using  the  assumed  KiVjOg)  of  0.01,  one  can  obtain 
an  approximation  of  the  act.lVjOg)  at  which  reaction  [11  should  commence  by  thermodynamic 
calculation  using  Gibbs  energy  data  provided  by  Yokokawa  et  al  (18),  viz., 

AG%ooc  =  IA.GV.  -  =  -RTInK  [2] 

(act.YVO^)*  X  (act.  ZrO,)  1 

K  =  _ 1  =  1.803x10’’  =  _ 

(act.VjOe)  X  (act.YjOg)  (4.5  x  10^)  x  (act.VjOg) 

or,  therefore  act.(V206)  1 .2  x  10®  for  reaction  Ill  to  occur 

This  value  of  act.(V206)  is  exceeded  by  the  act.(V206)  of  ~  10"*  found  for  NaVOj  at  850X  (see  below), 
and  therefore  it  is  consistent  that  YSZ  reacts  readily  with  molten  NaVO,  (Na20.V206),  but  not  with  the 
more  basic  NagVO^  (3Na20.V206),  as  found  in  the  next  section. 


Lewis  Acid/Base  Reaction  of  Ceramic  Oxides  with  Sodium  Vanadates 


Fig.  2  shows  the  reaction  behavior  of  a 
number  of  oxides  of  increasing  acidity  as  these 
are  heated  at  700-900®C  with  vanadium 
compounds  of  increasing  acidity;  i.e.,  NajVO^ 
(BNajO.VjOe),  NaVOj  (Na20.V206),  and  pure 
V2O6  (19).  Two  important  aspects  of  Fig.  2 
should  be  noted.  First,  the  reactions  are 
primarily  controlled  by  the  Lewis  acid/base 
nature  of  the  oxide,  with  no  reaction  occurring 
when  the  Lewis  basicities  are  approximately 
equal.  A  compound  such  as  NaVOg  can  react  as 
an  acid,  providing  VjOg,  with  an  oxide  of  the 
basicity  of  Y2O3,  but  as  a  base,  providing  NajO, 
with  oxides  of  the  acidity  of  Ge02  or  Ta206. 

Second,  the  ceramic  oxides  react  with 
either  the  Na20  or  V2O6  component  of  the 
sodium  vanadate  compounds,  and  not  with  the 
compounds  themselves.  The  role  of  the  sodium 
vanadate  compound  is  only  to  make  NajO  or 
V2O6  available  at  some  certain  level  of  activity. 
The  reaction  products  consist  of  multiples  of 
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Fig.  2  Lewis  acid/base  reaction  of  selected  oxides 
with  sodium  vanadates  and  V2O6. 


Na20  or  VjOg  with  the  various  oxides,  as  for  example,  2Na20.9Ge02,  Na20.2Ta206,  Y2O3.V2O6  (i.e., 
2  YVO4),  etc.  An  exception  to  this  rule  was  found  in  later  experiments  with  SC2O3,  which  reacts  with 
NaV03  to  form  3NaV03.Sc203  at  temperatures  below  880X  (20).  The  oxidation  state  of  the  different 
cations  does  not  change,  except  for  the  case  of  Ce02  where  Ce'**  goes  to  Ce®*  during  the  formation 
of  CeVO^  (see  below). 
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Previous  Studies  of  Vanadate-Sulfate  Thermochemistry 


Although  numerous  publications  on  vanadate-sulfate  hot  corrosion  can  be  found  in  the 
literature,  only  four  pertain  to  the  actual  thermochemistry  of  vanadate-sulfate  deposits.  These  are 
listed  in  Table  1 . 


TABLE  1 

Past  Studies  of  the  Thermochemistry 
of  Vanadate-Sulfate  Melts 

Tooic/Conclusions 


Thermodynamic  calculation  of  Na,V,S  deposition  In  gas 
turbines:  found  k(V206)  of  0.18  and  0.08  at  750*  and  900*C 
for  Na/V  =  2  and  8.7  x  10"*  and  3.6  x  10-^  pISOa),  resp. 

Electrochemical  study  of  NajO-VjOg  melts,  with  NajO  activity 
by  electrode  EMF,  VaOg  activity  by  Gibbs-Duhem  Integration; 
melts  found  to  be  highly  nonideal 

Calculation  of  NajVOg,  NaVOj  and  VjOg  concentrations  In  melts; 
vanadate-sulfate  melts  assumed  to  be  Ideal;  predicted  ~20% 
in  Na-V  melts  at  900®C  and  -log  Na20  of  15.3  —  18.7 

TG  A  study  of  SO3  equilibrium  with  NaVOg  at  800*C,  and  reaction 
of  melt  with  CeOj;  for  lowest  V2O6  levels,  found  k(V206)  ~  0.01 

Question  of  in  Vanadate-Sulfate  Melts 

In  the  older  boiler  corrosion  research  (1),  it  was  postulated  by  many  that  vanadium-rich 
deposits  contained  \/**/\/^*  couples  which  transported  oxygen  by  reactions  of  the  type: 

Na20.V204.5V206  +  Na20.6V206  (3) 

where  one  could  have 

Na20.V204.5V206  +  SO3  -*  Na20.6V206  +  SOj  (4) 

which  might  then  be  followed  by 

NajO.eVjOg  +  Fe -*  Na20.V204.5V206  +  FeO  (5) 

One  solution  proposed  for  vanadate-induced  corrosion  by  heavy  oils  in  boilers  was  the  use  of 
"low  excess  air",  or  LEA.  In  pilot  boiler  tests  (25)  using  low-quality  oil  containing  2.57  wt%  S,  53  ppm 
Na,  and  350  ppm  V,  the  vanadium  in  boiler-tube  deposits  forming  at  1300®F  (704*0  was  91%  V2O6, 
9%  V2O3  and/or  V2O4  with  15%  excess  air,  but  only  24%  V2O6  with  76%  V2O3  and/or  V2O4  at  1% 
excess  air.  In  the  latter  instance,  there  was  significantly  lower  deposit  formation  and  boiler-tube 
corrosion,  presumably  because  of  the  higher  melting  points  and  physical  nature  of  V2O3  and  V2O4. 
These  data  therefore  indicate  that  VjOj  or  V2O4  can  be  found  in  boiler  deposits.  However,  this  occurs 
under  nearly  reducing  combustion  conditions,  and  also  where  there  is  the  possibility  of  V2O6  reduction 
by  tube  metal  oxidation  (the  1 300°F  tube  deposits  in  fact  consisted  of  about  75%  corrosion  product). 


Year  Authors/Reference 
1982  Luthra  &  Spacil,  Ref.  21 

1984  Mittal  &  Elliott,  Ref.  22 

1 989  Hwang  &  Rapp,  Ref.  23 

1 995  Reidy  &  Jones,  Ref.  24 
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Whether  V**  would  occur  in  significant  amounts  in  vanadate-sulfate  melts  on  TBCs  in  gas 
turbines  operating  at  high  air/fuel  ratios,  or  in  laboratory  studies  of  vanadate-sulfate  melts  under 
ambient  Oj  partial  pressures,  is  another  question.  In  these  cases,  the  conditions  are  oxidizing,  and  the 
ceramic  oxides  involved  (ZrOj,  YjO,,  etc.)  are  stable  to  very  low  plOj),  except  for  CeOj  which,  as  the 
pure  oxide  at  1000“C,  begins  to  be  reduced  to  CeO,.83  at  atm  of  Oj.  Calculations  by  Hwang 

and  Rapp  (23)  indicate  that  NajSO^-GO  mol%  NaVO,  melts  at  gOO^C  and  ambient  plOj)  should  contain 
-20%  W**  over  the  log  act.lNajO)  range  of  -15.3  to  -18.7.  However,  this  contradicts  the  findings 
of  Mittal  and  Elliott  (22)  who  conclude  that  vanadium  in  liquid  NajO-VjOg  exists  only  as  over  the 
measured  log  act.(Na20)  range  of  -13  to  -20  at  852*C  (see  below).  They  support  their  conclusion  by 
citing  Block-Bolton  and  Sadoway  (26)  who  found  no  in  liquid  VjOe  at  1200K  (927*C)  for  p(Oj)  > 
0.06  atm.  Also,  it  appears  that  NajO  stabilizes  V®*  in  VjOg,  as  proposed  by  Mittelstadt  and 
Schwerdtfeger  (27)  who  observed  that  whereas  NajO.SVjOg  begins  to  reduce  to  at  10'^  atm  p(02) 
at  1000®C,  NajO.VjOg  is  stable  down  to  10"®  atm  p(02). 

As  another  consideration,  when  Ce02  reacts  with  the  V2O6  component  of  vanadate-sulfate 
melts  (see  below),  the  reaction  goes  as, 

2  CeOj  (s)  -f-  VjOg  (I)  -  2  CeVO^  (s)  +  1/2  Oj  (g)  [6] 

A  reduction  occurs  in  reaction  16],  but  one  can  not  be  certain,  a  priori,  whether  it  is  the  cerium  cation 
(4*  to  3*)  or  the  vanadium  cation  (5*  to  4*)  that  is  reduced.  To  answer  this  question,  Reidy  and 
S wider  (28)  examined  CeV04  by  X-ray  absorption  spectroscopy  (XAS)  and  showed  that  cerium  is  in 
the  Ce®^  state  in  CeVO^,  and  vanadium  thus  in  the  V®*  state.  Taken  overall,  the  literature  data 
indicate  that,  for  the  oxidizing  conditions  postulated,  is  unlikely  to  be  significant  in  the 
thermochemistry  of  vanadate-sulfate  melts  relating  to  TBC  ceramic  hot  corrosion. 

Thermodynamic  Modeling  of  Acid-Base  Melts 

The  activity  of  a  component  of  a  solution  or  melt  is  defined  according  to 

act.  =  K  X  X  [7] 

where  act.  is  the  component  activity,  y  is  the  activity  coefficient,  and  X  is  the  mol-fraction  of  the 
component.  For  ideal  solutions,  the  activity  coefficient  is  unity,  and  the  activity  of  a  component  is 
simply  equal  to  its  mol-fraction.  Almost  all  melt  solutions  are  nonideal,  however,  and  in  that  case,  y 
changes  with  the  overall  composition  of  the  solution.  To  describe  how  y  varies  with  melt  composition, 
solution  scientists  often  adopt  "solution  models"  which,  following  theoretical  or  empirical  reasoning, 
are  intended  to  predict  the  dependence  of  y  on  melt  composition. 

Many  solution  models  have  been  developed  since,  in  different  types  of  melts,  the  melt 
components  interact  in  different  ways,  which  must  be  taken  into  account.  For  oxide  systems  which 
involve  strong  acid-base  chemical  interactions  (such  as  at  present),  the  best  model  appears  to  be  the 
Ideal  Mixing  of  Complex  Components  (IMCC)  model  developed  by  Bonnell  and  Hastie  (29).  This  model 
is  capable  of  describing  the  very  large  changes  in  the  ratio  of  formal  concentration  to  thermodynamic 
activity  (10’®  or  more)  that  can  occur  in  strong  acid-base  reactions,  as  demonstrated  by  Bonnell  and 
Hastie  for  complex  slags  containing  such  components  as  Na20,  Si02,  K2O,  and  AI2O3. 

The  IMCC  model  uses  Gibbs-energy-minimization  by  computer  to  determine  the  "free"  mol- 
fraction,  as  opposed  to  the  nominal  mol-fraction,  of  a  component  such  as  K2O,  when  that  component 
exists  in  the  melt  in  competing  equilibria  with  other  components  by  reactions  such  as, 

K2O  +  SiOj  K2Si03  [8] 
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KjSiOj  +  SiOj  ^  KjSijOg 

[9] 

KjSijOg  +  2  SiOj  fi  KjSuOg 

1101 

It  is  thus  akin  to  the  NASA-developed  computer  calculation  of  the  equilibria  of  molecular  species  in  high 
temperature  gases,  and  also  to  the  calculations  made  for  the  NajO-VaOs-SOa  system  by  Luthra  and 
Spacil  (21).  The  IMCC  model  assumes  ideal  mixing  of  the  various  "product  molecules"  once  they  have 
been  formed  by  the  chemical  reactions  indicated;  that  is,  there  is  assumed  to  be  no  other  type  of 
interaction  between  the  product  molecules,  after  they  are  formed,  that  would  yield  significant 
additional  Gibbs  energy  of  solution  for  the  melt. 

This  raises  the  question  of  whether  such  "product  molecules"  actually  exist  in  melts.  For 
gases,  techniques  such  as  high  temperature  mass  spectrometry,  infrared  spectroscopy  or  gas 
chromatography  can  be  used  to  identify  gas  phase  molecules,  and,  if  conditions  are  such  that  kinetics 
do  not  interfere,  the  thermodynamically-predicted  gaseous  species  generally  are  found,  and  In  the 
approximate  quantities  calculated.  With  melts,  however,  it  is  not  possible  to  make  an  equivalent 
molecular  identification  of  melt  species.  About  the  best  information  available  is  obtained  with  binary 
metal  melts  such  as  Sb-Mg,  where  plotting  the  excess  solution  stability  (d^GVdx,*,  where  G^  Is  the 
excess  integral  molar  free  energy  of  a  solution)  against  melt  composition  may  reveal  peaks  at  certain 
stoichiometries  which  indicate  the  formation  of  melt  complexes  such  as  MggSbj  (30).  Bonnell  and 
Hastie  (29)  approach  this  problem  by  considering  that  "—  the  liquid  components  are  not  necessarily 
independent  molecular  or  ionic  species,  but  serve  to  represent  the  local  associative  order."  The 
goodness  of  fit  between  solution  model  prediction  and  experimental  results  obtained  by  Bonnell  and 
Hastie  appears  to  justify  their  hypothesis,  although  it  remains  difficult  to  envision  the  arrangement  of 
atoms  in  the  melt  that  this  would  imply. 

For  the  purposes  of  the  present  Report,  the  discussion  in  this  section  has  mostly  only 
philosophical  implications,  since  no  solution  model  for  the  NajO-VjOg-SOs  system  is  proposed  from  our 
NRL  research.  The  NRL  results  do  contain,  however,  substantial  experimental  data  (see  below)  that 
could  be  used  in  the  development  of  such  a 
solution  model. 

Study  of  Na,0-V,0c  Melts  bv  Mittal  and  Elliott 

Mittal  and  Elliott  (22)  consider  sodium 
vanadate  deposits  to  be  mixtures  of  vanadium 
oxide  (VjOg)  and  sodium  oxide  (NajO)  produced 
by  oxidation  of  Na,V-impurities  in  fuel  during 
combustion.  Our  NRL  research  has  taken  the 
same  approach,  but  with  the  effect  of  engine 
gas  SO3  added  in.  Since  the  Mittal  and  Elliott 
paper  in  some  ways  represents  the  "starting 
point"  for  our  NRL  work,  their  paper  will  be 
discussed  in  detail  here. 

Mittal  and  Elliott  made  electrochemical 
measurements  at  757-937”C  of  the  activity  of 
NajO  in  melts  having  a  VjOg  mol-fraction  range 
of  0.5  to  0.9;  that  is,  totally  molten  mixes  from 
NajO.VjOg  (i.e.,  NaVOj)  up  to  Na20.9V206  in  the 
VjOg-rich  half  of  the  NajO-VjOg  phase  diagram 
(Fig.  3).  These  composition  and  temperature 


Fig.  3  Phase  diagram  for  V^Og-NajCOs  (NajO). 
Ref.  31. 
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ranges  are  generally  thought  to  be  among  the 
most  corrosive  in  terms  of  vanadate  hot 
corrosion.  The  corresponding  activities  for  VjOg 
were  calculated  by  integration  of  the  Gibbs- 
Duhem  equation. 

The  activities  of  NajO  and  VjOg  found  by 
Mittal  and  Elliott  for  molten  NajO-VjOg  mixes  at 
852‘’C  (1 1 25K)  are  plotted  in  Fig.  4.  At  0.5  mol- 
fraction  of  VjOg  (that  is,  for  NajO.VjOg  or 
NaVOj),  the  log  of  NajO  activity  is  -13.13,  and 
that  of  VjOg  is  -3.62,  giving  a  log  activity 
product  of  -16.75.  This  value  can  be  compared 
with  thermodynamic  calculation  for  the  reaction, 

2  NaVOa  (I)  NajO  (I)  -t-  VjOe  (I)  11 1] 

where  the  log  activity  product  of  Na20  and  VjOg 
is  found  to  be  -1 6.93,  using  thermodynamic  data 
by  Yokokawa  et  al  (18)  and  Barin  (32).  This  is 
good  agreement  which  supports  the  Mittal  and 
Elliott  results  as  being  correct. 


Fig.  4  Log  act.lNajO)  and  log  act.lVjOe)  in  the 
NajO-VjOe  binary  system  at  1 1 25  K  (852*0.  From 
Ref.  22. 


The  data  of  Mittal  and  Elliott  are  plotted 

in  Fig.  5  with  NajO  and  VjOg  as  the  melt  components.  This  plotting  shows  VjOe  to  have  a  strong 
negative  deviation  from  ideal  solution  behavior.  The  V^Og  activity  coefficient  for  XlVjOg)  =  0.5  (i.e., 
in  molten  NaVOj)  is  ~4.8  x  10"^,  as  obtained  by  dividing  the  VjOg  activity,  log  -3.62,  by  the  VjOg  mol- 
fraction  concentration,  0.5.  The  Mittal  and  Elliott  data  can  also  be  recalculated  and  plotted  to  reflect 
NaVOj  and  VjOg  as  being  the  melt  components.  This  yields  Fig.  6,  where  VjOg  continues  to  exhibit 
a  negative  deviation  from  ideal  solution  behavior.  Fig.  6  indicates  therefore  that  NaVO,  and  VjOg  do 
not  mix  in  a  totally  ideal  fashion,  as  assumed  in  the  calculations  by  Hwang  and  Rapp  (23).  As  a  final 


Na^O.V^Og  mol-froction  V^Og 

Fig.  5  Solution  behavior  of  NajO-VjOg  at 
852°C,  taking  Na^O  and  V^Og  as  the  melt 
components. 


NaVOj  mol— fraction  VjOg  '^2^5 


Fig.  6  Solution  behavior  of  NajO-VjOg  at 
852°C,  with  NaVOj  and  VjOg  as  the  melt 
components. 
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side  point,  notice  that  there  are  no  obvious  inflections  in  the  experimental  log  act.lNajO)  curve  of 
Mittal  and  Elliott  (Rg.  4)  such  as  might  reflect  any  influence  in  the  melt  behavior  resulting  from  the 
compounds,  NajO.SVjOg  (NaVjO,)  and  NaaO.eVjOg  (NajVijOj,),  which  occur  within  the  NajO-VjOe 
composition  range  studied  (Rg.  3). 

Thermooravimetric  Equilibrium  of  SO,  with  NaVO.^  Melts 

The  mutual  interaction  between  NajO, 

VgOg  and  SO,,  as  occurs  in  engine  vanadate- 
sulfate  deposits,  was  investigated  by  a  thermo- 
gravimetric  technique  in  which  SO,  at  increasing 
partial  pressures  was  equilibrated  with  molten 
NaVO,  at  800*C.  The  reaction  is  described  by, 

2  NaVO,  (!)  +  SO,  (g) 

VjOe  (I)  +  Na2S04  (I)  (121 

Reaction  [1 21  was  identified  by  Luthra  and  Spacil 
(21)  as  being  the  major  reaction  in  determining 
the  composition  of  vanadate-sulfate  engine 
deposits  under  marine  gas  turbine  conditions. 

Also,  we  have  shown  earlier  that  whereas  CeO, 
does  not  react  chemically  with  pure  molten 
NaVO,,  an  overpressure  of  SO,  yields  VjOg  by 
reaction  (121,  with  CeV04  then  forming  by 
reaction  between  CeO,  and  VjOg  (19),  Note 
that  one  starts  with  pure  NaVO,  in  reaction  (121 
and  then  essentially  adds,  by  using  initially  very  ^  Arrangement  of  specimen  melt  planchet  (FT) 
low  partial  pressures  of  SO,,  small  Increments  of  within  furnace. 

VjOg  and  Na2S04  to  the  original  NaVO,,  The 

thermochemical  behavior  of  the  resultant  melts  would  therefore  be  expected  to  show  a  smooth 
transition  from  the  thermochemical  properties  found  by  Mittal  and  Luthra  for  0,5  mol-fraction  NajO- 
VjOg  (NaVO,),  at  least  for  the  initial  low  VjOg-  and  Na2S04-content  mixes. 

Our  thermogravimetric  analysis  (TGA)  furnace  system  used  a  two-stage  gas  dilution  arrange¬ 
ment  which  employed  electronic  mass  flow  gas  controllers  to  give  SO2  partial  pressures  of  down  to 
1 0'*  bar  In  the  furnace  input  air  (33).  The  SOj/air  input  mixture  was  equilibrated  over  Pt  at  temperature 
in  the  furnace  (Fig.  7),  with  the  S0,-S02-air  mixture  then  passing  over  the  NaVO,  melt  (50  mg)  which 
was  contained  in  a  shallow,  balance-suspended  Pt  planchet.  After  an  initial  weight  loss  of  ~2%  (the 
NaVO,  contained  about  2%  water),  the  system  showed  negligible  weight  change  (<  0.2  mg)  for 
periods  in  excess  of  50  hrs  under  slowly  flowing  (50  ml/min)  air  dried  over  Drierite™. 

Attainment  of  SO3  equilibrium.  As  the  0(50,)  was  lowered,  the  time  required  to  reach 
equilibrium  became  progressively  longer,  so  that  for  p(SO,)  ^  1 0'*,  1 00  hrs  or  more  would  be  needed 
for  equilibrium  by  simple  exposure  to  the  test  p(SO,).  Our  procedure  therefore  was  to  bring  the  NaVO, 
to  both  above  and  below  equilibrium  weight  (in  separate  measurements)  at  a  higher  piSO,),  and  then 
return  the  system  to  the  test  p(SO,)  and  monitor  the  wgt  gain  vs.  time  slopes  (positive  or  negative), 
which  were  back-extrapolated  to  identify  the  true  equilibrium  weight  (33),  The  time  for  equilibrium 
also  was  found  to  be  proportional  to  the  weight  of  NaVO„  taking  approximately  twice  as  long  for  100 
mg  as  50  mg  NaVO,.  This  experience  raises  the  question,  at  least  in  the  author's  mind,  of  how 
massive  amounts  of  melt  (up  to  25  gms)  in  a  narrow  electrochemical  cell  could  become  equilibrated 
with  low  pISOjI's  (^  10"*  atm)  in  a  reasonable  time.  Two  possible  answers  come  to  mind:  first,  many 
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of  the  electrochemical  experiments  were  at  900“C  rather  than  800‘’C  as  here,  and  the  100®C  higher 
temperature  may  significantly  speed  up  SO3  diffusion;  second,  our  melt  is,  at  least  initially,  pure 
NaVOj,  whereas  most  electrochemical  experiments  involving  SO3  equilibrium  use  melts  consisting 
mostly  of  sulfates,  and  it  may  be  that  SO3,  perhaps  as  SjOy",  diffuses  more  rapidly  through  sulfate- 
based  melts  than  vanadate-based  melts. 

Thermodynamic  Analysis  of  SO,-NaVO,  Eouilibrium  Data 


The  equilibrium  weight  gain  data  for  50 
mg  (0.41  mmols)  of  dried  NaV03  under  different 
pCSOj)  at  800®C  are  presented  in  Fig.  8.  The 
weight  gain  behavior  was  reversible  and  repro¬ 
ducible  (to  within  ±0.2  mg),  which  confirms 
that  SO3  is  absorbed,  or  given  off,  without  any 
permanent  change  in  the  NaVOg  melt.  Conver¬ 
sion  of  the  NaV03  completely  to  Na2S04  and 
VjOg  by  reaction  [12]  would  yield  an  equilibrium 
weight  gain  of  1 6.4  mg. 

From  the  weight  gain  and  stoichiometry 
of  reaction  [12],  one  can  calculate  the  mol- 
fraction  of  VjOg  experimentally  formed  at  each 
plSOg),  as  listed  in  Table  2.  And,  since  the 
activity  coefficient  equals  unity  for  ideal 
solutions  and  therefore  act.  =  X,  the  "ideal" 
mol-fraction  can  be  calculated  from  thermo¬ 
dynamic  data  following  the  procedure  outlined  in 
reaction  (2).  Using  the  thermodynamic  data 
provided  by  Luthra  and  Spacil  (21)  for  the 
reagents  and  products  in  reaction  [12],  one 
obtains. 


Fig.  8  Weight  gain  behavior  of  NaVOg  equilibrated 
at  800°C  with  the  indicated  SO3  partial  pressures. 


act. (V ,0k)  X  act.(Na7S0.) 

K  =  91.61166=  -  2  6»  <>20^41 

act.(NaV03)^  X  p(S03) 

which,  by  transposing  pISOa)  and  setting  act.  =  X,  reduces  to  a  quadratic  equation  in  the  unknown, 
X,  that  can  be  solved  by  numerical  approximation  on  the  computer  for  each  plSOg), 


x^ 

91.61166  xp(S03)  =  -  [14] 

(1-2x)^ 

The  "ideal"  mol-fractions  (or  "ideal"  activities)  for  VjOg  so  obtained  are  listed  in  Table  2.  Division  of 
the  ideal  VjOg  mol-fraction  by  the  experimental  VjOg  mol-fraction  then  yields  an  approximate  ylVjOs), 
listed  in  Table  2  as  "ylVjOgl-Approx.  1"  for  each  p(S03),  or  alternatively,  for  each  melt  composition. 

However,  approximation  1  for  KiVjOg)  is  not  correct  since  the  actual  value  for  act.lVjOg)  has 
not  been  determined.  Approximation  1  treats  NaVOs,  Na^SO^  and  VjOg  as  all  behaving  ideally,  which 
is  almost  certainly  not  true.  Luthra  and  Spacil  (21),  for  example,  considered  NaVOj  and  Na2S04  to  be 
ideal,  but  V^Og  to  be  nonideal  because  of  its  stronger  interactions  with  the  melt.  If  VjOg  is  assumed 
nonideal,  while  NaV03  and  Na2S04  are  ideal,  then  equation  [14]  can  be  modified  to  give,  using  the 
experimentally-determined  mol-fractions,  a  second  approximation  for  ylV^Og)  according  to. 
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HVjOgl-Approx.  2 


91 .61166  X  plSOg)  X  XCNaVOj)* 
X(NajS04)  X  X{V206) 


(151 


The  values  for  K(V206)-Approx.  2  are 
listed  in  Table  2,  with  both  approximations  being 
plotted  in  Rg.  9.  Inspection  of  Rg.  9  suggests 
that  KiVjOgl-l  is  probably  closer  to  correct  at  the 
higher  pISOj),  while  k(V206)-2  is  more  nearly 
correct  at  the  lower  pISO,)  values.  This  may  be 
partially  explained  along  the  following  lines. 

Note  in  equation  [13]  that  at  the  lowest  pCSOs), 

NaVOj  is  almost  pure,  and  its  activity  thus 
essentially  1 .  Under  these  circumstances,  K  x 
pISO,),  a  known  value,  becomes  equal  to  the 
activity  product,  i.e.,  act.(V206)  x  act.(Na2S04). 

Determining  act.(V206)  is  then  a  matter  of 
correctly  apportioning  the  activity  product.  In 
Approx.  1,  all  species  are  treated  as  ideal,  and 
act.(V206)  =  act.(Na2S04),  with  act.(V206)  being 
therefore  simply  1  /2  of  the  activity  product.  In 
Approx.  2,  where  Na2S04  (I)  but  not  V2O6  (D  is 
assumed  ideal,  act.(Na2S04)  =  X(Na2S04),  and 
division  of  the  activity  product  by  X(Na2S04) 
from  Table  2  yields  a  lower  activity  for  the 
"nonideal"  VjOg  Species.  In  actuality,  Na2S04  (I) 
is  probably  not  totally  ideal,  and  act.(V206) 

should  fall  somewhere  between  the  Approx.  1  and  Approx.  2  values.  One  might  expect  also  that,  as 
their  mol-fraction  concentration  increases  at  high  pfSOa),  V2O6  and  Na2S04  would  exhibit  activity 
coefficients  closer  to  unity  (cf.  Rg.  5),  making  Approx.  1  more  nearly  correct.  Therefore,  for  purposes 
of  a  later  comparison  (see  below),  an  averaged  k(V206)  was  calculated,  which  is  shown  as  ylN/jOgl-Avg 
in  Rg.  9.  The  calculated  average  is  weighed  toward  k(V206)-2  at  the  lower  pISOj),  and  kIVzObM  at 
the  higher  plSOg). 


p(SOj)  bar 

Fig.  9  Values  for  k(V206),  Approxs.  1  and  2,  and 
K(V206)-Avg,  plotted  against  plSOg). 


TABLE  2 

Experimental  Results  and  Calculated  Approximations 
for  KiVjOg)  in  NaVOj  Equilibrated  with  SO3  at  800°C 


plSOa) 

Wgt.  gain 

Exptl. 

"Ideal" 

K(V206)- 

K(V206)- 

K(V206)- 

in  bar 

in  mo 

XIVjQsl 

XlY2Qa 

Aonrox.  1 

Anorox.  2 

Avg. 

5.6  X  10  ® 

2.2 

0.0670 

0.00225 

0.0336 

8.6  X  10-^ 

0.0055 

4.1  X  10'^ 

2.8 

0.0853 

0.00603 

0.0707 

3.6  X  10  ® 

0.023 

3.3  X  10^ 

3.8 

0.116 

0.0167 

0.144 

1.3  X  10® 

0.069 

2.0  X  10-® 

5.8 

0.177 

0.0392 

0.221 

2.5  X  10  ® 

0.14 

1.0  X  10^ 

7.6 

0.232 

0.0799 

0.344 

4.9  X  10  ® 

0.26 

5.8  X  10"* 

11.0 

0.335 

0.157 

0.469 

5.1  X  10® 

0.41 

9.9  X  10-^ 

12.5 

0.381 

0.186 

0.488 

3.6  X  10  ® 

0.49 
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Note  that  the  k{V206)  values  in  Table  2  depend,  of  course,  upon  the  thermochemical  data 
selected.  If  data  incorporating  values  provided  more  recently  by  Yokokawa  et  al  (18)  are  used,  the 
value  of  K  is  raised  from  91 .61166  to  168.8  which,  for  example  at  plSOg)  =  5.6  x  10  ®  bar  changes 
KiVjOsl-Approx.  2  for  from  8.6  x  1 0'*  to  1 .6  x  1 0  *. 

Reaction  of  Candidate  Stabilizing  Oxides  with  SOo-NaVO. 

The  equilibrated  SOg-NaVOg  TGA  technique  was  used  to  investigate  the  reaction  of  MgO,  Y  O  , 
SC2O3  and  In^Oj  (33),  CrjO,  and  SnOj  (34),  and  CeOj  (24)  with  the  NajO-VjOg-SOa  melt  systern,  with 
the  mix  ratio  of  0.41  mmols  NaVOj  and  0.1  mmols  of  oxide  normally  being  employed.  As  summarized 
in  Rg.  lOa-lOg,  the  results  serve  to  rank  the  potential  vanadate-sulfate  hot  corrosion  resistance  of 
candidate  stabilizing  oxides  for  ZrOj.  They  also  provide  information  about  the  nature  of  the  different 
reactions  that  the  various  oxides  undergo  with  vanadate-sulfate  melts.  These  can  be  broken  down  into 
categories  as  described  below: 


Fig.  10  (a-g)  Weight  gain  vs.  plSO,)  at  800°C  for  NaVO,  (0.41  mmols)  plus  various  candidate  ZrO,- 
stabilizing  oxides  (0.10  mmols). 


^  ,n7  of  soluble  product.  This  is  exemplified  by  MgO  (Fig.  10a),  where  even  at  a  p(SO,) 

of  10  bar,  the  MgO-NaVOj  mix  shows  a  higher  weight  gain  than  for  pure  NaVOj,  even  though  the 
amount  of  NaVOj  is  the  same  in  each  case.  The  additional  weight  gain  must  come  from  the  uptake 
of  SO3  by  either  MgO  or  Na^O,  since  vanadyl  sulfate  is  ruled  out  at  such  low  p(S03).  Low  melting 
mixtures  of  VjOg-fZMgO.NaiO)  or  VjOs-lMgO.Na^O)  can  form  (35),  but  these  would  yield  no  change 
in  weight.  Also,  thermodynamic  calculations  indicate  that  the  displacement  of  Na^O  by  MgO  by,  e.g.. 
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MgO  +  NajO.VjOs  -*  NajO  +  MgO.VjOe,  is  highly  unfavored.  Solid  MgS04  also  appears  not  to  be 
possible,  since  a  critical  p(S03)  of  1 .6  x  1 0'*  bar  is  Indicated  for  MgO  (s)  +  SO3  -*  MgS04  (s)  at  800*0. 
However,  If  soluble  MgS04  is  formed,  then  appreciable  MgS04  could  be  produced  even  at  p(S03)  ~ 

1 0^  bar.  Moreover,  since  such  a  system  would  have  two  degrees  of  freedom  by  the  phase  rule,  some 
MgS04  (1)  can  be  formed  no  matter  how  low  the  piSO,);  i.e.,  there  will  be  no  one  critical  pISOj)  for 
onset  of  reaction.  The  hypothesis  of  a  soluble  MgS04  product  is  consistent  also  with  the  MgO-ZrOj 
TBC  burner  rig  tests  (12)  burning  Na,V,S-containing  fuel,  where  major  amounts  of  MgS04  but  only 
minor  amounts  of  Mg  vanadate  were  found  on  the  TBC  surface.  At  the  higher  pfSOj)  values  in  Fig. 
1 0a,  there  is  an  upturn  in  weight  gain  which  presumably  reflects  the  onset  of  magnesium  vanadate 
formation,  since  the  reaction,  2MgO  +  VjOs  =  2Mg0.V206  is  thermodynamically  predicted  to  be 
possible  at  800®C  at  act.CVjOg)  =  8x10"*  bar.  This  could  not  be  confirmed  unequivocally  by  x-ray 
diffraction,  however,  which  in  this  case  gave  complex  and  diffuse  patterns  (common  for  Na/Mg 
vanadate-sulfate  phases  from  melts).  Nonetheless,  taken  overall,  the  weight  gain  behavior  in  Fig.  10a 
indicates  that  melt-soluble  MgS04  is  formed  at  low  plSOg)  in  the  MgC-NaVOs-SOg  system,  but  that  Mg 
vanadate  may  be  produced  at  high  plSOg)  (i.e.,  at  high  VjOg  melt  activities). 

Metal  vanadate  formation  once  critical  V/)f  activity  reached.  This  behavior  is  shown  by  CeOj 
(Rg.  10b),  In203  (Fig.  10c),  and  ScjO,  (Fig.  lOd),  although  the  case  for  SC2O3  is  complicated  by  the 
occurrence  of  3NaV03.Sc203.  With  Ce02  and  In203,  the  weight  gain  curve  is  just  the  same  as  pure 
NaVOg  (in  other  words,  the  oxide  is  as  chemically  inert  as  the  Pt  planchet  itself)  up  to  a  critical  plSOg), 
where  the  melt  V2O6  activity  is  raised  (via  reaction  1 2)  sufficiently  that  the  oxide  begins  to  react  with 
V2O6  to  form  the  vanadate.  The  weight  gain  for  Ce02  and  In20s  is  then  almost  vertical,  reflecting  the 
insolubility  of  the  vanadate  product.  Since  the  rare  earth  oxides  and  their  vanadates  are  high-melting 
compounds  with  excellent  x-ray  diffraction  patterns,  one  can  show  by  x-ray,  as  for  Ce02  and  CeV04 
(see  below),  that  only  oxide  exists  below,  and  only  vanadate  above,  the  critical  pISOs),  or  actually  V2O6 
activity. 


With  SC2O3,  the  weight  gain  curve  for  Sc203-NaV03  Is  below  that  of  NaV03  up  to  the  critical 
pISOj).  Although  not  originally  understood  (32),  this  behavior  results  because,  at  temperatures  below 
880*C  (20),  SC2O3  reacts  with  NaVOg  to  give  the  weak  compound,  3NaV03.Sc203.  This  compound 
reduces  the  activity  of  Na20,  and  therefore  less  SO3  is  taken  up  than  with  pure  NaVOg  of  equivalent 
weight.  The  weight  gain  rise  leans  somewhat  from  the  vertical  in  this  case,  which  possibly  reflects 
the  "breaking"  of  the  3NaV03.Sc203  compound  as  piSOg)  increases.  Also,  our  Initial  X-ray  analysis 
was  in  error,  with  SCVO4  being  correctly  identified  as  the  only  scandium  compound  above  the  critical 
p(S03),  but  SC2O3  as  the  compound  below,  when  in  fact  it  was  3NaV03.Sc203.  The  difficulty  in 
distinguishing  between  SC2O3  mixed  with  NaVO,  and  3NaV0,.Sc203  arises  because  all  of  the  strong 
peaks  for  SC2O3  (x-ray  standard  pattern  JCPDS  42-1 463)  are  also  contained  in  the  3NaV03.Sc203  x-ray 
pattern  (JCPDS  30-1241). 

Direct  reaction  with  pure  NaVO^.  In  Fig.  lOe,  Y2O3  is  sufficiently  basic  that  it  reacts  directly 
with  molten  NaV03  to  produce  YVO4,  which  leaves  the  melt  enriched  in  NajO.  Therefore,  when  this 
system  is  equilibrated  with  even  very  low  p(S03),  there  is  a  high  weight  gain  because  of  the  pre¬ 
existing  "reservoir"  of  NSjO. 

Another  type  of  direct  reaction  occurs  with  Cr203  (Fig.  1 0f)  which  reacts  with  pure  NaV03  and 
reduces  the  activity  of  NajO,  probably  by  Na2Cr04  formation.  This  then  yields  less  uptake  of  SO3  than 
would  be  found  with  pure  NaV03  alone,  and  over  a  wide  piSOg)  range. 

Inert  to  the  NaVOySOj  system.  Finally,  Sn02  (Fig.  lOg)  is  chemically  inert  to  the  NaV03-S03 
system,  with  the  weight  gain  being  essentially  identical  to  that  for  NaV03  alone  up  to  at  least  1  x  1 0  ^ 
bar  of  SO3.  The  resistance  of  Sn02  to  reaction  with  V2O6  may  be  one  reason  for  its  effectiveness  in 
passivating  fluid  cracking  catalysts  against  deactivation  by  vanadium  impurities  in  crude  oil  (7). 
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Determining  Act.(V-.Ocl  without  Assumption  as  to  Act.(Na-SOJ 


In  calculating  act.lV^Og)  from  the  NaVOa-SOj  equilibrium  weight  gain  data  via  equation  (131, 
one  is  ultimately  dealing  with  the  activity  product.  acUV^O^)  x  act.(Na2S04),  and  an  assumption 
concerning  the  Na2S04  solution  behavior  (ideal  or  nonideal,  and  what  activity  coefficient?)  must  be 
made  to  obtain  act.{V206).  However,  an  alternative  derivation  of  act.(V206),  or  more  specifically,  of 
K(V206),  that  is  not  dependent  upon  act,(Na2S04)  can  be  made  on  the  basis  of  reaction  (61, 

2  CeOj  (s)  +  VjOe  (I)  2  CeV04  (s)  +  1/2  O2  (g)  (6] 

In  this  case,  Ce02  and  CeV04  are  pure 
solids  of  low  solubility  (24),  and  plOj)  is  the 
atmospheric  oxygen  pressure  of  0.21  bar  so,  at 
800*C,  the  activity  of  V206(l)  is  fixed.  Using  the 
thermodynamic  data  for  CeOjfs),  CeV04(s)  and 
V206(l)  provided  by  Yokokawa  et  al  (18),  the 
critical  V2O6  activity  necessary  for  reaction  (61 
can  be  calculated  to  be  6.2  x  10"^. 

The  problem  is  to  identify  the  VjOg 
concentration  at  which  reaction  (61  just  begins. 

This  can  be  done  by  examination  of  Fig.  11, 
which  shows  that  the  weight  gain  step-up 
resulting  from  reaction  (6)  commences  at  a 
pISOa)  of  3  X  10"®  bar,  where  the  melt  weight 
gain  is  3.0  mg,  corresponding  to  formation  of 
0.091  mol-fraction  of  VjOg  (and  Na2S04).  If  the 
theoretical  V2O6  activity  for  Reaction  (61  of  6.2 
X  1 0'*  is  divided  by  the  experimental  VjOg  mol- 
fraction  of  0.091  from  Rg.  11,  a  VjOg  activity 
coefficient  of  7  x  10®  results.  This  k(V206)  ^‘3-  11  Weight  gain  behavior  for  CeOj-NaVOs 

agrees  well  with  the  KlVjOgl-Avg.  of  ~0.01  defining  critical  plSOg)  for  Ce02 -*  CeV04  reaction, 
indicated  in  Fig.  9  for  pfSOj)  =  4  x  10'^  bar, 

even  though  HVzOg),  in  the  present  case,  is  determined  on  the  basis  of  a  different  experimental 
reaction,  and  without  regard  to  the  Na2S04  or  NaVOj  solution  behavior. 

As  is  well  known,  calculations  such  as  above  depend  critically  upon  the  accuracy  of  the 
thermodynamic  data.  Assuming  an  accumulated  ±  1 0  kJ  error  in  AG*  for  reaction  (61  causes  the 
derived  k(Y206)  to  vary  from  2.1  x  10*  to  2.2  x  10'*.  This  error  bar  is  shown  for  the  datum  point  of 
K(V20g)  ==  7  X  10®,  X(V206)  =  0.091  which  has  been  included  as  representative  of  the  present 
experiment  for  the  k(V206)  comparison  in  Fig.  14  below.  It  is  not  an  insignificant  error,  but  fits  well 
within  the  data  band.  Good  agreement  is  thus  seen  between  two  independent  sets  of  experiments 
(i.e.,  S03-NaV03  equilibrium,  and  Ce02  CeV04  reaction),  in  which  two  different  sets  of 
thermodynamic  data  (Luthra  &  Spacil,  and  Yokokawa  et  al)  were  used.  Error  in  the  thermodynamic 
data  is  therefore  unlikely  to  seriously  invalidate  the  present  results,  especially  in  their  role  of  serving 
to  indicate  the  general  range  of  V2O6  nonideality  in  vanadate-sulfate  melts. 

V,0g  Titration  using  CeVQ.  Formation  as  "Indicator*  to  Determine  v(V..O.;l 

Reaction  (6)  is  analogous  to  the  metallurgical  reaction, 

M  (s)  +  1/2  O2  (g)  MO  (s)  [16] 
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which  has  great  significance  in  metallurgy.  By  reaction  [161,  one  can,  for  a  given  temperature, 
determine  either  A,G*  for  MO,  or  the  critical  plOj)  for  metal  oxidation/reduction,  provided  that  the  other 
is  known.  Moreover,  reaction  [161  is  totally  independent  of  the  gaseous  environment.  Large  quant¬ 
ities  of  inert  gases  such  as  nitrogen  or  argon  can  be  present,  and  there  can  be  other  reactions  in  which 
O2  engages  (e.g.,  CO  +  1  /2  02  =  COj)  but,  if  obeyed  as  written,  reaction  [1 61  will  nonetheless  remain 
independent. 


For  the  case  of  reaction  [61,  a  critical 
VjOg  (I)  activity  of  6.2  x  10^  is  required  at 
800*C,  according  to  calculations  based  on  the 
thermodynamic  data  of  Yokokawa  et  al  (18),  for 
the  reaction  to  proceed, 

2  CeOj  (s)  +  VjOe  (I) 

2  CeVO^  (s)  +  1/2  O2  (g)  [61 

If  obeyed  as  written,  reaction  [61  is  independent 
of  the  melt  composition,  and  the  VjOg  (I)  activity 
required  for  reaction  [61  will  always,  at  800*0, 
be  6.2  X  10"^  regardless  of  melt  composition. 
However,  the  VjOg  activity  coefficient  can,  and 
most  probably  will,  change  with  melt 
composition.  Therefore,  by  finding  the  mol- 
fraction  of  V2O6  at  which  reaction  [61 
commences,  one  can  potentially  determine 
K(V206)  for  virtually  any  melt  system  in  which  an 
appropriate  V2O6  activity  for  reaction  [61  to 
proceed  can  be  attained. 


Mmols  of  VjOj  Added 


Fig.  12  Mmols  VjOg  required  for  CeV04  formation 
at  800®C  vs.  NaV03-Na2S04  melt  composition. 


This  idea  was  tested  using  melts 

consisting  of  1  mmol  of  Ce02  mixed  with  4  mmol  blends  of  vanadate-sulfate  of  four  different 
compositions,  viz.,  1 00%  NaVOg,  NaVOs-d  0  mol%)Na2S04,  NaV03-(25  mol%)Na2S04,  and  NaV03-(50 
mol%)Na2S04.  A  series  of  varying  amounts  of  were  added  to  the  melt  mixes,  which  were  then 
heated  at  800*C  for  2  hrs  in  quiescent  air  and  x-rayed  to  determine  the  amount  of  Ce\/04  formed.  The 
relative  amount  of  CeV04  produced  was  determined  via  the  expression. 


[CeV04(200)  Peak  Hgt.l 

%  CeV04  =  100  X  - f - - - 

[CeV04{200)  PeakHgt.1  +  [Ce02(111)  Peak  Hgt.l 
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which  is  an  approximation  of  the  integrated  intensity  ratio  formulation  commonly  used  to  determine 
the  relative  amounts  of  crystal  phases  by  x-ray  diffraction  (36).  By  extrapolation  of  the  percent  of 
CeV04  formation  to  zero  (Fig.  13),  one  can  ascertain  the  mmols  of  V2O6  that  must  be  added  to  raise 
the  melt  VjOe  activity  to  the  point  at  which  CeV04  formation  begins,  and  Ce02(s),  CeV04(s),  VjOed), 
and  atmospheric  O2  are  all  in  equilibrium;  that  is,  when  reaction  [61  is  obeyed. 

The  results  yielded  from  Fig.  12  are  tabulated  in  Table  3,  where  the  calculated  k(V206)  is  6  x 
1 0^  for  NaVOj-d  0  mol%)Na2S04  to  which  0.1 031  XlVjOg)  had  been  added.  This  XlVjOg)  is  near  the 
0.0853  X(V206)  obtained  by  equilibrating  NaVOj  under  4.1  x  lO"'  bar  of  SO3  (cf.  Table  2),  and  where 
K(V206)-Approx.2  was  3.6  x  10  ®.  To  obtain  a  better  comparison,  the  data  in  Table  2  were  linearly 
interpolated  to  0.1031  X(V206),  which  yielded  a  value  for  K(V206)-Approx.2  of  9.1  x  10  ®.  This  latter 
j/'IVjOs)  of  9.1  X  10  ®  agrees  well  the  k(Y206)  of  6  x  10  ®  from  Table  3,  especially  considering  that  a 
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linear  interpolation  was  applied  to  the  Table  2  data.  A  datum  point  of  HVjOg)  =  6  x  10  *.  XfVjOg)  = 
0.1031  with  an  error  bar  corresponding  to  ±10  kJ  error  in  which  is  still  the  largest  likely 

error,  is  included  in  Fig.  14  as  representative  of  this  experiment.  As  shown  in  the  drawn  square  in  Rg. 
1 4,  the  three  independent  experimental  techniques  all  give  essentially  the  same  KiVjOg)  for  VjOg  in  the 
0.1  mol-fraction  concentration  range. 


TABLE  3 

Millimoles  of  VjOg  Required  for  CeV04  Formation 
as  a  Function  of  Na2S04  Content  of  the  Melt  at  BOOX 


Nominal  melt 
composition 

100%  NaVOs 
NaVOj-IIO  mol%)Na2S04 
NaV03-(25  mol%)Na2S04 
NaV03-(50  mol%)Na2S04 


Mmols  V2O6 
for  CeV04 

Mol-fracti 

formation 

of  V.O. 

0.51 

0.1131 

0.46 

0.1031 

0.38 

0.0868 

0.29 

0.0676 

Final 

mol-fraction  Act.  Coeff.’ 
of  Na,SO,  of  V,0. 

0.0  0.0055 

0.0897  0.0060 

0.2283  0.0071 

0.4662  0.0092 


’  Calculated  for  act.(V206)  =  6.2  x  10"^ 


When  V2O6  is  added  to  melts  containing 
significant  Na2S04,  the  back-reaction  of  Reaction 
11 2]  can  occur,  with  NaV03  produced  in  the  melt 
and  SO3  evolved.  In  a  flowing  air  stream,  where  ^ 
SO3  is  continually  removed,  there  can  be  a  major  S 
extent  of  back-reaction,  with  50  mol%  Na2S04- 
NaVOa  being  converted  ~  80%  to  NaV03  in  24  « 

hrs  at  750®C  II 9).  In  quiescent  furnace  air,  the  ° 
degree  of  back-reaction  is  much  less.  Also  the  ^ 
reaction  kinetics  are  slow  in  comparison  to  most  -t 
chemical  reactions.  In  the  present  case,  the  o 
CeOj  is  finely  ground  together  and  thoroughly 
mixed  with  the  NaV03-Na2S04-V206  blend  before  > 

firing.  Upon  coming  to  temperature,  there  is 
rapid  reaction  to  form  stable  CeV04,  and  since 
the  Ce02  is  in  excess  of  VjOg,  no  residual  V2O6 
is  left  to  undergo  reaction  with  Na2S04  in  the 
melt  (which  is  heated  only  2  hrs).  Error  because 
of  Na2S04  reaction  with  residual  V2O6  should 
therefore  be  negligible. 


Fig.  13  Activity  coefficient  of  VjOb  vs.  mol- 
When  the  VjOg  activity  coefficients  from  fraction  of  Na2S04  at  800”C. 

Table  3  are  plotted  against  the  mol-fraction  of 

Na2S04,  an  essentially  straight  line  is  obtained  (Fig.  1 3).  Theoretical  reasons  as  to  why  the  line  should 
be  straight  have  not  been  developed.  However,  Figs.  12  and  13  contain  two  points  which  have  not 
been  explicitly  shown  before,  and  which  are  perhaps  not  intuitive.  First,  increasing  the  Na2S04  content 
of  the  melt  causes  k(V206)  to  increase;  that  is,  V2O6  reaction  with  Ce02  (and  presumably  other  ceramic 
oxides)  occurs  with  less  V2O6  required  in  the  melt  (or  vanadium  in  the  fuel).  Second,  small  differences 
in  the  V2O6  activity  coefficients  are  not  negligible.  A  change  in  ylVjOg)  of  from  5.5  x  10  ®  to  9.2  x  ■®, 
although  at  first  sight  insignificant,  in  fact  corresponds  to  a  nearly  2X  decrease  from  0.51  mmols  to 
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0.29  mmols  in  the  amount  of  VjOg  required  to  reach  the  melt  VjOg  activity  at  which  CeOj  will  react 
or  "corrode".  In  the  service  world,  this  could  mean  having  to  buy  a  fuel  of  half  the  vanadium  content 
at  considerably  more  expense. 

SUMMARY  AND  CONCLUSIONS 


Comparison  with  V,0c  Activity  Coefficients  from  the  Literature 


The  known  VjOg  activity  coefficients 
that  have  been  reported  for  sodium  vanadate, 
and  sodium  vanadate-sulfate  mixtures,  are 
shown  in  Fig.  14.  The  comparison  is  only 
qualitative,  since  the  activity  coefficients  were 
measured  with  different  melt  conditions  and 
temperatures.  The  Mittal  and  Elliott  data  point 
is  calculated  from  their  measurement  of  the 
activity  of  VjOg  in  0.5  mol-fraction  NajO.VjOg 
(i.e.,  NaVOg)  at  SBO'C.  The  two  Luthra  and 
Spacil  data  points  were  reported  for  melts 
having  a  NaA/  ratio  of  2,  and  equilibrated  under 
pISOjI's  of  approximately  9  x  10"^  and  4x10"^ 
atm  at  750*C  and  900'’C,  respectively.  Although 
Hwang  and  Rapp  (23)  originally  postulated 
sodium  sulfate-vanadate  solutions  to  be  ideal,  a 
later  paper  by  Rapp  and  Zhang  (37)  describes 
the  equilibration  of  Na2S04-V206  melts  of  varying 
Na/V  ratio  under  a  fixed  piSOj)  of  2.2  x  1 0'*  bar 
at  900®C,  and  concludes  that  the  melts  are  "-- 
close  to  an  ideal  salt  solution  although  a  small 
deviation  seems  to  be  present."  No  activity 
coefficient  for  VjOe  is  given  in  their  paper,  but  in 
a  review  of  an  NRL  manuscript  on  vanadate- 


Mol-fraction  VjOj 

Fig.  14  Comparison  of  VjOg  activity  coefficients 
from  various  works  as  plotted  vs.  mol-fraction  of 
melt  VjOe- 


sulfates,  a  reviewer  who  seemed  quite  knowledgeable  about  the  work  stated  that  a  VjOg  activity 


coefficient  of  0.79  was  indicated  by  the  Rapp  and  Zhang  results.  So,  despite  its  unusual  derivation. 


this  KlVjOg)  is  plotted  in  Rg.  14  for  completeness  of  the  comparison. 


Although  somewhat  scattered,  the  data  in  Fig.  14  exhibit  a  general  trend  of  increase  in  kIV^Ob) 
as  the  mol-fraction  of  VjOg  increases.  Despite  the  difference  in  temperature,  there  is  an  essentially 
smooth  transition  between  the  KlVjOg)  of  ~5  x  10"*  (or  10'**)  indicated  by  the  Mittal  and  Elliott  data 
for  pure  NaVOj  at  850®C,  and  the  kIV^Ob)  of  ~5  x  10'*  determined  from  NRL  results  for  an  800°C 
NaVOj  melt  containing  0.067  mol-fraction  each  of  NajSO^  and  VjOb.  This  is  consistent  with  our 
expectations,  considering  that  at  this  point  only  small  increments  of  VjOb  have  been  added  to  pure 
NaVOj.  There  is  also  reasonably  good  agreement  with  the  Luthra  and  Spacil  VjOb  activity  coefficients 
of  0.08  and  0.1 8  in  the  middle  regions  of  the  plot,  and  with  the  0.79  VjOb  activity  coefficient  of  Rapp 
and  Zhang  for  Na/V  melts  under  high  SO3  partial  pressures. 


The  three  datum  points  within  the  drawn  square  are  of  particular  importance.  As  described 
above,  these  points  were  determined  by  three  independent  experimental  techniques,  viz.,  simple 
equilibrium  of  NaVOs  under  SO3;  the  reaction  CeO^  +  VjOb  -*  CeV04  +  1/2  O2  where  the  V2OB  melt 
activity  was  increased  by  raising  pISO,);  and  the  same  reaction  where  V2O6  was  simply  "titrated"  into 
the  melt.  Two  different  sets  of  thermodynamic  data  were  also  employed.  Good  agreement  is  seen 


between  the  three  values  for  ^(VjOg)  so  obtained,  which  gives  confidence  that  our  experiments  and 
data  are  essentially  correct. 

Other  considerations 

Effect  of  Na/V ratio  in  melt  or  engine  deposit  In  addition  to  pISOj),  the  other  principal  variable 
in  the  Na20-V206-S03  system  is  the  NaA/  ratio.  Information  on  the  effect  of  the  HaN  ratio  is  contained 
in  our  VjOg  "titration"  experiments  where  blends  of  different  NaVOj-NajSO^  ratios  were  investigated. 
For  example,  NaVOs-ISO  mol%)Na2S04  has  a  NaA/  ratio  of  3,  and  could  be  considered  as  approxi¬ 
mating  a  melt  that  was  originally  Na20-(50  mol%)NaV03  (or  Na3V04),  but  in  which  the  NajO 
component  has  been  sulfated.  Such  an  assumption  is  not  unreasonable,  since  NajO  has  a  considerably 
higher  affinity  for  SO3  than  for  V2O6,  as  evidenced  by  the  fact  that  the  800*C  NajSO^  dissociation 
activity  product,  act.(Na20)  x  act.(S03),  is  9.8  x  10'*\  whereas  the  eguivalent  activity  product  for 
NaV03,  act.(Na20)  x  act.lVjOg),  is  1.7  x  10'^®.  Taking  this  model  to  be  indicative  even  if  not  totally 
correct,  one  can  interpret  the  data  in  Fig.  13  as  being  for  melts  of  NaA/  ratios  of  1,  1.22,  1.67  and 
3.0,  respectively,  wherein  HVjOg)  is  found  to  progressively  increase  from  5.5  x  10’®  to  9.2  x  10'®,  or 
by  67%.  Compared  to  the  large  changes  in  ylVjOe)  produced  by  increasing  pCSOj),  this  is  relatively 
small,  which  suggests  that  KiVjOg)  is  not  highly  dependent  upon  the  melt  Na/V  ratio,  at  least  over  the 
concentration  ranges  covered. 

Vanadate-sulfate  melts,  idea!  or  nonideal?  The  data  in  Fig.  14  demonstrate  that  vanadate- 
sulfate  melts  are  nonideal,  and  in  fact  follow  a  very  common  solution  behavior  wherein  the  melt 
components,  when  present  in  high  concentration,  behave  as  "nearly  ideal",  but  when  in  dilute 
concentration,  exhibit  large  negative  deviations  from  Ideal  behavior,  and  have  small  activity 
coefficients. 
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